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Summary: In recent experiments the non-equilibrium distribution func- 
tion f{E, U) in diffusive Cu and Au quantum wires at a transport voltage 
U shows scaling behavior, f{E, U) — f{E/eU), indicating a non- Fermi liq- 
uid interaction with non-vanishing T — scattering rate. The two-channel 
Kondo (2CK) effect, possibly produced by degenerate two-level systems, is 
known to exhibit such behavior. Generalizing the auxiliary boson method 
to non-equilibrium, we calculate f{E, U) in the presence of 2CK impuri- 
ties. We show that the 2CK equations reproduce the scaling form f{E/eU). 
For all measured samples the theoretical, scaled distribution functions co- 
incide quantitatively with the experimental results, the impurity concen- 
tration being the only adjustable parameter. This provides a microscopic 
explanation for the experiments and, considering that no other mechanism 
producing the scaling form is known to date, lends strong evidence for the 
presence of degenerate two-level defects in these systems. The relevance of 
these results for the problem of dephasing in mesoscopic wires is discussed. 



1 Introduction 

Two-level systems (TLS) have been known to exist in disordered solids for a long 
time. Their signatures have been observed in the anomalous thermodynamic 
properties of glasses, which may be explained by a flat distribution of level split- 
tings. In metals, slow two-level fluctuators are evidenced by telegraph noise in the 
conductance A new physical phenomenon, the two-channel Kondo (2CK) ef- 
fect H, ^, arises when a fast, energetically (nearly) degenerate TLS is embedded 
in a metal, where the local impurity degree of freedom couples to the contin- 
uum of conduction electrons via an exchange interaction. The channel degree of 
freedom, conserved by this interaction, is the magnetic conduction electron spin, 
which is always degenerate in the presence of time reversal symmetry. 
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The 2CK effect exhibits striking non-Fermi hquid behavior |Q, ^, ^ with a 
non- vanishing single-particle scattering rate at the Fermi energy ep for tem- 
perature T = 0, logarithmic behavior of the linear specific heat coefficient and 
the susceptibility, and a non-analytic, universal correction to the electronic den- 
sity of states, AN{E)/{B^/T) — h{x), where _B is a non-universal constant, 
X := \E — sf\/T and h{x) is a universal scaling function with h[x) = ^/x for 
a; 1 [0, 1^. Signatures of the latter have been observed as zero bias conduc- 
tance anomalies (ZBA) of nano point contacts [|[ |l^ |ll|. As suggested in Ref. 
p2| , the finite single-particle scattering rate at the Fermi energy would provide a 
natural explanation for the saturation of the dephasing time r^p at low tempera- 
tures, which has recently been observed in magnetotransport measurements 
of weak localization in disordered wires. 

However, the actual existence of the 2CK effect in nature has remained a 
controversial issue, partly because the ZBA could also be explained qualita- 
tively 0, ahhough not quantitatively by the Al'tshuler-Aronov 
diffusion-enhanced Coulomb interaction ||17|, partly because the physical realiza- 
tion of degenerate TLS is poorly understood. Therefore, it is essential to develop 
unambiguous methods for detecting 2CK physics in mesoscopic systems. 

Recently, it has been demonstrated in a landmark experiment performed by 
the Saclay "Quantronics" group that unique information about the elec- 
tronic interactions in a metal may be extracted from the distribution function 
f{E^ U) of quasiparticles with energy E when the system is driven far away from 
equilibrium by a transport voltage U: The shape of /(i?, U) is determined by 
the energy dependence of relaxation processes which tend to equilibrate the sys- 
tem. It has been found that in diffusive Cu and Au nano-wires the distribution 
function obeys a scaling form, /(i?, U) — f{E/eU), implying a non-Fermi liquid 
interaction |l§ . 

In this article we discuss the close relation between this peculiar scaling prop- 
erty, which we will call "-E/et/ scaling" for brevity, and non- Fermi liquid behav- 
ior. It is then shown that the 2CK effect obeys E/ell scaling and reproduces 
quantitatively the measured distribution functions, while any other type of in- 
teraction is ruled out. This provides the strongest case to date for the physical 
realization of the 2CK effect in Cu and Au nano-wires induced by TLS. We 
briefly discuss its relation |Q to the problem of the dephasing time saturation 
at low temperatures . 

2 Experiment and non-Fermi liquid signature 

The distribution function f{E,U) was measured []l8| in an experimental setup 
where a non-equilibrium current was driven through a Cu nano- wire contacted by 
two reservoirs at chemical potentials and eU, respectively (Fig. |l|). In addition. 
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a superconducting Al tunneling junction was attached at a position x along the 
wire, the Al slab being in equilibrium with itself. For a voltage V across the 
junction the tunneling current is given by 



Jt'U 



J dE[f{E) - r{E + eV)]NcuAE) Nsc{E + eV) , (2.1) 



where t is the (energy independent) tunneling matrix element, f°{E) is the Fermi 
distribution function in the superconductor, and Ncu, Nsc denote the density of 
states in the wire and in the superconductor, respectively. Since for the voltages 
used in the experiment Ncu is flat and the BCS density of states N^c is measured 
independently, the non-equilibrium distribution f{E) in the Cu wire can be 
extracted from this expression. The electronic transport in the wire is diffusive 
with diffusion coefficient D. Length L and thickness d of the wire are such that 
the Fermi surface is three dimensional and the Coulomb and phonon scattering 
times are large compared to the electronic diffusion time tjj — L^/D through 
the wire, so that equilibration due to these processes can be neglected jl^. In 
this situation, assuming purely elastic scattering, one expects the distribution 
function at a given position x along the wire to be a linear superposition of the 
Fermi functions in the reservoirs, 

feiast^j^^ [/) = (1 - ^)r{E + eU) + y°{E) , (2.2) 
since there is no energy exchange within the electron system or between the 



electron system and the lattice. Eq. (2.2) is a solution of the diffusive Boltz- 



mann equation ( |3.6D [ |2l[ , when the collision integral vanishes (see section 3). 
This situation is to be distinguished from the hot electron regime, where local 
equilibration occurs due to inelastic processes p2[ . 

The measured distribution functions showed rounding of the Fermi steps as 
compared to Eq. ( |2.2| ) and obeyed scale invariance with respect to the transport 
voltage [/, f{E,U) — f{E/eU), when U exceeded a certain low energy scale, 
eU - 0.1 meV Deviations from scaling were observed again for voltages 
larger than a high energy scale E^ — O.bmeV. The latter may be explained by 
reservoir heating effects or by the electrons coupling to additional degrees of 
freedom at high energies. We now deduce the non-Fermi liquid signature from 
the scaling property. The latter implies that the equation of motion for f{E), the 
Boltzmann equation and, as a consequence, the inelastic single-particle collision 
rate 1/t{E) are scale invariant. Assuming a (yet to be determined) two-particle 
potential V{e) with energy transfer e, 1/r is given in 2nd order perturbation 
theory as Q 

1 1 ,„,Q f , f , ,,~ , -/ E 



t(E) T{E/eU) 



NcuiOf j j de'\Vie)?E[l^,±j,^) . (2.3) 
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Figure 1 Experimental setup for measuring tlie non-equilibrium distribution 
function in a metallic nano-wire. 



Here is a combination of distribution functions / guaranteeing that there is 
only scattering from an occupied into an unoccupied state. Therefore, the ex- 
perimental results about the scaling property of / imply that F depends only on 
the dimensionless energies as displayed in Eq. (2.3). Demanding scale invariance 
for 1/r with respect to eC/, i.e. making the frequency integrals dimensionless, 
implies a characteristic energy dependence of the interaction and the single par- 
ticle scattering rate, V{e) cx 1/e 18| and 1/t{E) cx —\n{E/Eo), for energies 



less than Eq. These infrared singularities indicate a breakdown of Fermi liquid 
theory within the 2nd order perturbation theory argument applied here. 



3 Two-channel Kondo effect and scaling in non-equilibrium 

The strong infrared divergence of the two-particle potential, V{e) oc 1/e, de- 
duced above, is not explained by any conventional interaction, including the 
Al'tshuler-Aronov interaction It must, therefore, be generated by an infi- 
nite resummation of logarithmic terms obtained in perturbation theory due to 
the presence of a Fermi edge. In this section the 2CK effect is briefly reviewed. 
We then show that the effective electron-electron vertex, mediated by a 2CK 
impurity, has a 1/e divergence and calculate the resulting distribution functions 
away from equilibrium. 

The 2CK effect arises whenever a local, energetically degenerate two-level 
degree of freedom (pseudospin r ~ ±1/2) is coupled to a system of two identical 
conduction electron bands or channels via a pseudospin exchange interaction, 
which, however, conserves the channel degree of freedom. As for the physical 
realization of the pseudospin, it has been suggested M that it might arise from 
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Figure 2 Schematic snapshot of the Kondo screening cloud. The bold arrows 
represent the impurity (pseudo)spin. a) Single-channel Kondo effect. Conduction 
electrons and local impurity spin form a collective singlet ground state with en- 
tropy ^(r = 0) = 0. b) Two-channel Kondo effect. Each of the two identical 
conduction electron bands form a screening cloud (black and grey), so that the 
net pseudospin of the combined clouds is not 0; pseudospin singlet formation is 
frustrated. The ground state is not unique, resulting in a non-vanishing zero-point 
entropy ^(r = 0) = fcBln\/2. 

the degenerate positions of an atom in a symmetrical double well potential, 
from a rotational degree of freedom of a lattice defect or group of atoms, or 
from sliding kinks on screw dislocations in a lattice iQ. The channel degree 
of freedom is then the magnetic electron spin a = ±1/2, which is necessarily 
degenerate because of time reversal symmetry (Kramers degeneracy). 

In order to describe non-equilibrium situations, it is useful to represent this 
system in terms of the low-energy physics of a SU(2)pscudospinXSU(2)channci An- 
derson impurity model in the Kondo limit, i.e. in terms of a doubly degenerate 
local level Ed < 0, coupled via a hybridization v to two identical conduction 
channels, with conduction electron operators c^.^^., cj.^^. The dynamics of the 
two-level degree of freedom (pseudospin) is described by a fermionic field opera- 
tor dj.^ , with the restriction that the local level must not be doubly occupied at 
any time. The latter is implemented exactly by decomposing the local operator 
into auxiliary fermion and boson bg. operators, dl^. = f^bg., supplemented by 
the operator constraint = Er flfr+T.a ^l^s = 1 H- Thus, we describe the 
2CK system by the Anderson hamiltonian in auxiliary particle representation, 

H = Y1 +EdJ2 frfr + « T.(frb,C,r^ + h.C.) . (3.4) 

k(TT T kar 

Note that the boson field bg transforms with respect to the adjoint representa- 
tion of SU(2)channoi (denoted by the index ct). By integrating out bg in the Kondo 
regime ((Er/r/r) — 1^ l-^l l-^'dl), H is reduced to a 2CK model where the 
exchange coupHng J = [■up/ji^dl is mediated by the boson propagator and, as 
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Figure 3 Leading contributions to the effective electron-electron vertex induced 
by a 2CK impurity. Solid, dashed, and wiggly lines correspond to conduction 
electron, pseudofermion, and slave boson propagators, respectively. 



required, conserves the channel degree of freedom a. The above auxihary parti- 
cle representation is known to describe the infrared behavior of the 2CK effect 
correctly ||2^ already on the level of a leading-order selfconsistent perturbation 



theory in the hybridization v, the non-crossing approximation (NCA) |27 . It 



also allows for the application of standard Green's function methods and, in 
particular, is readily generalized to non-equilibrium by means of the Keldysh 
technique. 

Physically, the presence of two identical conduction electron channels leads 
to an overscreening of the impurity pseudospin below the Kondo temperature 
Tk = (2r/7r)exp(-7r|£;rf|/2r), T = 7r|t>piVc„(0), as elucidated in Fig. |. This 
results in a non-vanishing zero-point entropy S{T = 0) = A:BlnV2 Q and is 
closely related to the finite single-particle decay rate 1/t{E = 0, T = 0) > 
predicted by the 2CK effect. 

Below the Kondo scale Tk, the auxiliary particle propagators Gf, Gb are 
known to exhibit infrared power-law behavior in equilibrium, Gf{Lu) oc to^^f , 

Gb{i^j) oc uj^"'' for lo^Tk where the exact [26: relation af + ah — 1 



is characteristic for the multichannel Kondo effect (a/ = M/{M -\- N), ab 
N/ (M -|- iV) for pseudospin degeneracy N and channel number M) . Using power 
counting arguments, it follows that each fermion (boson) propagator contributes 
a power —a/ (— ttb) and each frequency integral contributes a power 1 to the 
frequency dependence of the effective electron-electron vertex 7(0;), shown in 
Fig. ||. Thus, the 2CK electron-electron vertex depends on the frequency transfer 
u as 7(w) oc ijj'^~'^(°'f+°'>') — ijj-i ^j^g perturbative level, this implies scale 

invariance of the inelastic scattering rate 1/t{E) with respect to the transport 
voltage U, as discussed in section 2. A self-consistent resummation of these 
terms in non-equilibrium shows that the scale invariance persists up to single- 
particle energies E of the order of eU . It may be shown that the range of elJ 
values for which scaling behavior is obeyed is bounded from below by the Kondo 
temperature Tk- For large bias voltage, deviations from scaling occur when eU 
becomes comparable to the band cutoff or to the local level \Ed\. 
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Figure 4 Scaled non-equilibrium distribution functions for a short (L = 1.5/xm, 
diffusion const. D = 65cm^/s) and a long [L — 5.0jj,m, diffusion const. D = 
45cm Vs) Cu nano-wire. Black curves: Experimental data for f/ in a range be- 
tween O.OSmV and 0.3mV in steps of O.OSmV, taken from Ref. [Q. The two curves 
showing deviations from scaling are at small voltages, 0.05mV and O.lmV, respec- 
tively. Light curves: Theory in the scaling regime, eU S> Tk- The 2CK impurity 
concentration dmp obtained from the fit is approximately 8 • 10~®/(lattice unit 
cell). 
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E/eU 

Figure 5 Deviations from scaling. Experimental data are tlie same as in Fig. 4, 
upper panel. Open circles represent the fit of the theory to the first experimen- 
tal distribution function which shows deviations from scaling at small voltage, 
eU=0.1meV. Since the 2CK impurity density Cimp has been determined as the 
only adjustable parameter in the scaling regime (Fig. 4), for given transport volt- 
age U, the Kondo temperature Tk (or the temperature T, whichever is larger), 
is the only parameter that controls deviations from scaling at small voltage. It 
may be detetermined from the fit as Tk — IK. 

The remaining task is to calculate the non-equilibrium distribution function 
fx{E, U) at an arbitrary position x along the wire (compare Fig. |^) in the pres- 
ence of 2CK defects. In the regime where the inelastic mean free path imei, 
approximately given by the average spacing a between 2CK impurities, is small 
compared to the wire length, but large compared to the elastic mean free path 
£, L ^ £inei ~ a ^ the dynamics of a given 2CK impurity is determined 
by interactions with electrons probing many different spatial regions in the 
wire. It is, therefore, appropriate to first calculate the averaged quasiparticle 
distribution {E, U) as the solution of the impurity ensemble averaged Boltz- 
mann equation, subject to the boundary conditions fx=o{E, U) = f°{E -\- eU), 
fx=L{E^ U) = f°{E). The equation of motion for / is obtained by summing over 
all local quasiparticle momenta p and exploiting the fact that the current in the 
disordered system is diffusive. 
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Ue) = V -ME, U) = -D\/^p,{E) - -D\/, V UE, U). (3.5) 
^ — ' m — ' 

p p 

Here, denotes the component of the gradient in the direction along the wire. 
In the stationary case, the resulting diffusive Boltzmann equation reads |21 



VlUE.U) = --C2CK[ME,U),C,mp], (3.6) 

where C2CK is the collision integral due to 2CK scattering. For small impurity 
concentration Cimp is proportional to Cimp- Since the diffusion constant D is 
measured experimentally, Cimp is the only adjustable parameter of the theory in 
the scaling regime. 



In the experiment, the tunneling current Eq. (2.1) measures the local (un- 
averaged) distribution function of electrons in the vicinity of the junction. It 
is determined by the local stationarity condition that it must be equal to the 
distribution function of the 2CK impurity states. The results are shown in Fig. ^ 
and display scaling behavior for eU 3> Tk- It is seen that there is excellent quan- 
titative agreement between theory and experiment for all samples. Considering 
the fact that no other interaction producing E/eU scaling is known to date, 
this provides strong evidence for the presence of 2CK impurities in evaporated 
Cu nano-wires. From the fit of the theory to that experimental curve which for 
the first time shows deviations from scaling as the voltage is decreased, i.e. for 
elJ = O.lmeV (Fig. 5), the corresponding experimental low energy scale may be 
identified with the low energy scale of the theory, the Kondo temperature Tk- 
Thus we have Tk — IK, which is in rough agreement with earlier experimental 
results on TLS in Cu point contacts j|, |l^ . 



4 Relation to dephasing 

In equilibrium, the non- vanishing 2CK quasiparticle scattering rate 1/t{E) cross- 
es over to a pure dephasing rate as the quasiparticle energy approaches the 
Fermi surface, i?,T — > 0, since at the Fermi energy no energy exchange is pos- 
sible. One might, therefore, conjecture that degenerate TLS could be the origin 
of the dephasing time saturation observed in magnetotransport measurements 
of weak localization. This assumption is indeed supported by several coinci- 
dences between the dephasing time measurements and the results on the non- 
equilibrium distribution function; (1) The dephasing time extracted from 
magnetotransport experiments [ pO| , |l3| is strongly material, sample, and prepa- 
ration dependent. This suggests a non-universal dephasing mechanism, like TLS, 
which is not inherent to the electron gas. (2) The dephasing time in Au wires is 
generically shorter than in Cu wires . This is consistent with the fact that the 
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estimates for the TLS concentration Cimp, obtained from the fit of the present 
theory to the experimental distribution functions, is much higher in Au than 
in Cu wires ]30[| . (3) In Ag wires one observes neither dcphasing saturation nor 
E/eU scahng of the distribution function. This is consistent with the assumption 
that there are no 2CK defects present in the Ag samples |^ . 

5 Concluding remarks 

We have calculated the quasiparticle distribution as a function of the excita- 
tion energy E in diffusive nano-wires in the presence of 2CK impurities, when 
the system is driven far away from equilibrium by a finite transport voltage 
U. The present theory reproduces the experimental finding that in Cu and 
Au wires the nonequilibrium distribution function displays the scaling prop- 
erty f{E, U) = f{E/eU) for eU above a low-energy scale Tq. Within the theory. 
To is given by the Kondo temperature Tk or by T, whichever is larger. The 2CK 
impurity density is the only adjustable parameter of the theory in the scaling 
regime. The quantitative agreement between the present theory and experiment 
and the fact that the experimental scaling property is not even qualitatively 
explained by any other type of interaction, provide strong evidence for the ex- 
istence of 2CK impurities in Cu and Au wires, while there seem to be no such 
defects present in Ag wires. We have not attempted here to give a microscopic 
model for the physical realization of 2CK defects. For that purpose, it should 
be useful to perform numerical simulations [p4[ of dislocations in the respective 
materials. 
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